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Atomistic density functional theory (DFT) calculations of the capacitance between a metallic 
cylindric gate and a carbon nanotube (CNT) are reported. Results stressing the predominant effect 
of quantum capacitance in limiting or even enhancing screening properties of the CNT are shown. 
Other contributions to the quantum capacitance beyond the electronic density of state (DOS) are 
pointed out. Negative values of the quantum capacitance are obtained for low-density systems, 
which correspondingly over-screen the gate field. This unconventional behavior of the quantum 
capacitance is related to the predominance of the exchange contribution in the total electronic 
energy of the CNT. 



Theoretical investigation of the physical properties of 
realistic quasi one-dimensional electronic systems like 
CNTs becomes more and more attractive due to the 
promising perspectives they offer to nanotecnology. Es- 
pecially interesting are the electronic properties of a CNT 
when it is used as a quasi-one dimensional channel in 
nanoscaled field-effect devices 0, 0, 0- The reduced 
dimensionality of these systems offers the possibility of 
directly investigating the fingerprints of the many-body 
nature of the electron-electron interactions in the screen- 
ing phenomena. In this Letter we report on atomistic 
DFT calculations of the screening properties of a coax- 
ially gated CNT (10,0). The signature of many-body 
physics manifests itself at low carrier densities where the 
CNT (10,0) exhibits a singular screening behavior, turn- 
ing from a poorly screening electron gas to a classical 
metal and eventually resulting into an over-screening sit- 
uation. 

It is known that for a system of low-dimensional con- 
ductors the classical electrostatic concept of capacitance 
can no longer be applied. Luryi reported that a two- 
dimensional electron gas (2DEG) used as grounded mid- 
dle plate in a three plates capacitor is not able to screen 
the field emanating from the charges on one of the other 
plates; on the contrary, the field partially penetrates 
through the 2DEG and induces charge on the opposite 
external plate. This non-classical field penetration is due 
to the reduced screening properties of the 2DEG, and it 
changes the capacitance from a geometrical quantity to 
an electrochemical one 0. In a circuital approach, this 
means that new quantum capacitances will appear be- 
sides the geometrical ones to describe the non-classical 
(quantum) screening properties of such low-dimensional 
electronic systems |4(. 

Quantum capacitance is however a more general con- 
cept which is strictly connected to that of thermody- 
namic compressibility. If we refer to an interacting one- 
dimensional electronic gas, the compressibility K is de- 



fined as la K~ 



d Etot/dn 2 = n 2 dp/dn, where E t , 



jj,, and n are, respectively, the total ground-state energy 
(per unit length), the chemical potential and the elec- 
tronic density of the system. The quantum capacitance 
per unit length Cq is defined as 1/Cq = (\j e 2 )d\ij dn 0. 
Clearly the capacitance results to be proportional to the 
compressibility. Combining the expressions for Cq and 
for K one also finds L/Cq = l/e 2 p (ep) ■ Kq/K, where 
Po(cf) and Kq are, respectively, the DOS at the Fermi 
energy ep and the compressibility of the free electron 
gas. This result shows that neglecting electron-electron 
interaction (Kq/K = 1), the quantum capacitance can 
be simply obtained from the DOS at the Fermi level. An 
exact evaluation of quantum capacitance is however of 
fundamental importance for nanoelectronic devices such 
as CNT field-effect transistors (CNTFET). Gate modula- 
tion depends on the combined effect of the insulator and 
quantum capacitance. Up to now, however, calculations 
on CNTFETs have only considered the proportionality 
between Cq and the DOS neglecting electron-electron 
interactions effects, which may induce strong variations 
of Cq. Model Hamiltonian approaches show that, due 
to the exchange interactions, the compressibility of quasi 
one- and two-dimensional electronic systems may also be 
negative 0- A negative sign of the compressibility has 
also been found in capacitive measurements on 2DEG re- 
alized in GaAs quantum wells [ljj. The purpose of this 
work is to shed some light on the quantum capacitance 
issue for CNTs, clarifying its dependence on electron- 
electron interactions. We demonstrate that the quan- 
tum capacitance of a CNT is not just determined by the 
DOS. In other words, the quasi one- dimensional nature 
of a CNT, reflected in a small DOS, may not only result 
in a poor screening behavior. On the contrary, when the 
carrier density is sufficiently low, due to the exchange 
interactions the CNT over-screens the external gate field 
and its quantum capacitance assumes negative values. 

In Fig. ^the schematic section of the simulated coax- 
ially gated CNT (10,0) is shown. The self-consistent 
screened potential needed for capacitance calculations is 
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FIG. 1: Coaxially gated (10,0) semiconducting CNT. The 
bare and screened electrostatic potential are shown as a 
dashed and solid line respectively. Plots refer to the cross- 
section taken in the middle of the gate contact. Potential 
inside the CNT is clearly over-screened. 

obtained solving in real space a three dimensional Poisson 
equation. The metallic gate contact is treated as an ap- 
propriate boundary condition for the electrostatic prob- 
lem. The gate is surrounded by an insulator layer with 
dielectric constant e r = 3.9, which prevents the CNT to 
touch the metal in realistic field-effect devices. Charge in- 
jection inside the CNT is obtained using two semi-infinite 
portions of CNT (10,0) which act as source and drain 
reservoirs Dirichlet boundary conditions have been 
used on the Poisson calculation box correspondingly to 
the source and drain contacts; this introduces two metal- 
lic surfaces in the system (see Fig. UJ. 

Calculations have been performed using the non- 
equilibrium Green's function [l^ density functional 
tight-binding (gDFTB HI El) 

approach, which is based 
on a local density approximation (LDA) of the exchange- 
correlation (XC) functional. The XC contribution to the 
electronic energy is modeled as an on-site Hubbard en- 
ergy It is important to note that the value of this 
Hubbard parameter takes into account the attraction in- 
duced by the exchange contribution to the on-site energy, 
and results to be smaller than an on-site parameter which 
just corresponds to the Hartree repulsion energy. Calcu- 
lations have been performed using only p z orbitals for 



carbon atoms. This allow us to simulate a 21.3 nm long 
CNT (10,0) (corresponding to 50 unit cells). We stress 
that we obtain the same over-screening results simulat- 
ing a shorter CNT modeled with a complete sp 3 basis. 
To obtain density dependent results, we artificially dope 
the semiconducting CNT (10,0) varying the number of 
valence electrons per carbon atom. In this doping model 
any added extra charge is uniformly distributed over the 
whole CNT, leaving the band structure unaffected 

Several theoretical methods have been proposed in the 
literature to extract the compressibility of a low dimen- 
sional electron gas. Standard approaches consist on nu- 
merical evaluations of the second derivative of the total 
energy with respect to the density 0, 0] . Other meth- 
ods based on the evaluation of a static structure factor 
which takes into account finite width effects in two and 
one-dimensional electron gases have also been proposed 
0. In this work we use a capacitive model (see Fig. 
^1 to directly extract the sign and the magnitude of the 
quantum capacitance at an atomistic scale. A different 
atomistic approach based on the evaluation of a properly 
defined inductance matrix .5] has been also presented 

ua 

The circuital model shown in Fig. ^extends the equiv- 
alent circuit used in Ref. |4| to a one-dimensional system 
with a locally applied gate field. The series connection 
of the geometrical capacitances C* s and C£ s de- 
scribes the capacitive induction between the gate and 
those parts of the source (drain) metallic surfaces inter- 
nal to the CNT. Using two capacitances we can account 
for the fact that this capacitive induction is not direct, 
but it is mediated by the CNT. In other words, part of 
the charge that the gate field would induce on the inner 
parts of the source(drain) metallic surfaces is induced on 
the CNT. The quantum capacitance Cq just accounts 
for the capability of the CNT of localizing this charge 
on its surface. When there are no mobile electrons able 
to screen the gate field, we have Cq — and charge is 
induced only on the source(drain) metallic surfaces. On 
the other hand, in the classical limit (\Cq\ — ► oo) the 
CNT completely screens the gate field, hence no charge 
is induced on the source (drain) metallic surfaces inside 
the CNT itself. Finite values of Cq describe partial in- 
duction of charge on the tube. The capacitance C DSiext 
accounts for the charge induced by the gate field on those 
parts of the metallic source(drain) surfaces external to 
the CNT. If the length of the CNT is large with re- 
spect to both the gate length and the insulator thick- 
ness, the charge induced on the metallic source(drain) 
surfaces is small with respect to the charge induced on 
the CNT. In this situation we can neglect both C^ s int 
and C DS cxt . This approximation apply to our system, 
where C^ s int /\CQ\ < 10~ 4 and C r>S-ext is at least two 
orders of magnitude smaller than the series of CL, , 
and Cq in the explored range of carrier densities |18| . 
Starting from the simplified circuital model, we obtain 
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FIG. 2: Computed quantum capacitance per unit length of 
the CNT (10,0) as a function of both the Fermi level and the 
carrier density. Only the first subband of the CNT is filled. 
Solid line refers to DFT calculations. Dashed line is obtained 
just considering the DOS contribution to Cq. 

the following meaningful expression for Cq 

where 8Q is the charge induced on the CNT by a small 
applied gate bias SVq, n is the CNT carrier density and 
the carrier density value which causes the CNT 
to totally screen (TS) the gate field. Clearly a negative 
quantum capacitance is obtained at those carrier densi- 
ties for which the CNT is able to localize more charge 
than what strictly needed to completely screen the gate 
field. We note that C x DSint is just the geometric insu- 
lator capacitance of the cylindric capacitor. Neverthe- 
less, due to the strong fringing field effects it can not 
be analytically evaluated. We operatively compute it as 
5Q (n TS ) /SVq, referring to the limit situation in which 
the CNT completely screens the gate field and a classical 
cylindric capacitor is recovered. 

In Fig. El we report the cross-section of the screened 
and bare potentials taken in the middle of the gate con- 
tact. Both bare and screened potentials correspond to 
a gate bias 5Vg — 10~ 4 V and a bias of V applied to 
both the source and drain metallic contacts. The CNT 
(10,0) is lightly n-doped, resulting in a small carrier den- 
sity n = 4.7 • lO- 2 !" 1 . In this situation the CNT (10,0) 
over-screens the gate field, as clearly shown by the pres- 
ence of an over-screened negative electrostatic potential 
inside the CNT itself. We stress that the over- screening 
result does not depend on the particular doping model we 
have used. In fact, a negative quantum capacitance has 
been also obtained for an intrinsic metallic CNT (5,0) 
|l9|. The only advantage in using a n-doped semicon- 
ducting CNT is that we can continuously vary the carrier 
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FIG. 3: Electron-electron interaction contribution to the 
quantum capacitance of a CNT (10,0). Solid line refers to 
the DFT atomistic calculation. Dashed line is an analytical 
result for the ratio Cq dos /Cq hf obtained for a parabolic- 
confined quantum wire matching the CNT. 

densities inside the first conduction subband. 

The presence of an additional contribution to the CNT 
quantum capacitance beyond the DOS can be appreci- 
ated in Fig. Here we plot the quantum capacitance 
computed for several carrier density values which cor- 
respond to a given Fermi level position inside the first 
CNT subband. Solid line refers to DFT calculations 
(Cq dft ), while the dashed line is just obtained plotting 
c Qdos = e2 p( e F), where p(e F ) is the DOS of the CNT 
computed at the Fermi level. C'q dos is the limit which is 
recovered by the DFT calculation for large enough car- 
rier densities. As the density lowers, Cq dft tends to 
an infinite positive value. This singular behavior corre- 
sponds to the critical density n TS for which the CNT 
completely screens the external field. Below this criti- 
cal density the CNT accumulates more charge than what 
strictly needed for total screening, and an over-screening 
situation is generated. We note that a negative value of 
the quantum capacitance does not mean necessarily in- 
stability. The presence of a positive insulator capacitance 
with a smaller absolute value than Cq gives a positive 
value to the total capacitance, thus stabilizing the whole 
system 

The presence of a negative quantum capacitance is 
related to the attractive exchange energy contribution, 
which at low carrier densities becomes predominant pro- 
ducing a local minimum in the total energy and a change 
in the sign of its second derivative (i.e. in the sign of 
the compressibility K) |(|. This theoretical prediction is 
confirmed by results we obtain eliminating the XC en- 
ergy contribution in the Hubbard on-site parameter: in 
this case our calculations result in a positive screened po- 
tential inside the CNT, on which in fact a smaller charge 
than SQ (n TS ) is induced. 

As explained earlier, the ratio Cq dos / Cq dft is rep- 
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resentative of the ratio Kq/ 'Kdft, where Kr>pT is the 
DFT compressibility of the CNT. Analytical evaluations 
of the effects that electron-electron interactions have on 
the compressibility of various quantum wire models have 
been presented in the literature |9|,|20| ■ We observe that if 
the radius of the CNT we deal with is not too big respect 
to the orbital extension of the carbon atom, the CNT can 
be suitably approximated by a parabolic confined cylin- 
drical quantum wire. A negative compressibility regime 
has been observed in such systems in the Hartree-Fock 
(HF) approximation, where electron-electron interactions 
are included at the lowest order as an exchange energy. 
Using these analytical results we have evaluated the ratio 
Cq D os/ c Qhf, obtaining 



Cn„„ 
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2m*E F 
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2m* E F 



(2) 



The function / is the static form factor which takes 
into account width effects associated to the quasi one- 
dimensional nature of the electron gas and b is the width 
parameter of the quantum wire. The quantities Ep and 
to* are respectively the Fermi level and the effective mass 
of the electronic system, while a* — ft 2 /m*e 2 is the ef- 
fective Bohr radius defined with m*. Results obtained 
using a width parameter b = 4A and an effective elec- 
tron mass m* ~ 0.083 m e (matching respectively the ef- 
fective width and the first subband effective mass of the 
CNT (10,0)) are shown in Fig. [31 together with the ratio 
Cq D os /Cq dft obtained within our approach. Compar- 
ing our atomistic DFT results with the analytical cal- 
culations, it is important to observe that in our model 
the exchange energy is only included in a local approxi- 
mation, and that the Hubbard approximation to the XC 
functional takes also into account correlation energy con- 
tributions. Moreover, our atomistic model includes the 
Hartree interaction which is completely absent in a uni- 
form electron gas model due to the presence of a neutral- 
izing positive background. Respect to purely HF models, 
the correlation contribution decreases the inverse com- 
pressibility, while the direct Coulomb contribution in- 
creases it [20|. Nevertheless, the good qualitative agree- 
ment of results presented in Fig. [3] shows that the fun- 
damental relationship between electron-electron interac- 
tions and over-screening in a CNT is correctly predicted 
by both DFT and HF calculations. 

Results shown both in Fig. |3 and in Fig. refer to 
a range of densities for which only the first subband is 
occupied. Adding further carriers, the second subband 
of the CNT starts to be filled and a more complex be- 
havior of the quantum capacitance appears (namely, we 
observe again a transition to negative values when the 
second subband starts to be filled). 

To summarize, atomistic DFT simulations have re- 



vealed that a CNT can over-screen an external field pro- 
duced by a gate contact. Negative values of the CNT 
quantum capacitance have been obtained in a low density 
regime inside the first conduction subband. This demon- 
strate that the DOS is not the only contribution to the 
capacitance of a gated CNT, and that electron-electron 
interactions effects can assume a fundamental role in de- 
termining the screening properties of a gated nanoscale 
device and consequently its field-effect dependence. 
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